Russian  Federal  Nuclear  Center 


* 


All-Russian  Scientific  Research  Institute  of 
Experimental  Physics 


APPROVED  BY 


Deputy  Scientific  Leader 
V.N  .Mokho  v 
#$1995  r. 


STUDYING  POSSIBILITY  OF  X  RADIATION  GENERATION  AT  LOW 
DENSE  PLASMA  COMPRESSION  USING  HIGH  -  POWER  EMG. 

* 

(Final  report  under  subcontract 
F  6170894W0757) 


19980303  139 


* 


MIC  QUALITY  m-SPECTED  & 


607200  Arzamas-16,  Nizhny  Novgorod  region.,  RF. 

1995 


gjffSTRIBUTION  STATE^MShTA 

Approved  lor  public  release; 
Distribution  Unlimited 


Form  Approved  OMB  No.  0704-0188 


REPORT  DOCUMENTATION  PAGE 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson 
Davis  Hidhwav.  Suite  1204,  Arlinqton,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503. 

1 .  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

1995  Final  Report 

4.  TITLE  AND  SUBTITLE 

Studying  Possibility  of  X  Radiation  Generation  at  Low  -  Dense  Plasma  Compression  Using  High 
Power  EMG 

5.  FUNDING  NUMBERS 

F61 70894 W0758 

6.  AUTHOR(S) 

Prof.  Vladislav  Mokhov 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Institute  of  Experimental  Physics 

Prospekt  Mira  37 

Arzamas-16  Russia 

Russia 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

N/A 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

EOARD 

PSC  802  BOX  14 

FPO  09499-0200 

10.  SPONSORING/MONITORING 

AGENCY  REPORT  NUMBER 

SPC  94-4076 

11.  SUPPLEMENTARY  NOTES 

12a.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 

12b.  DISTRIBUTION  CODE 

A 

1  1 3.  ABSTRACT  (Maximum  200  words) 

This  report  results  from  a  contract  tasking  Institute  of  Experimental  Physics  as  follows:  Consider  the  physical  scheme  of  an  arrangement  for 
x-ray  generation  by  quasi-spherical  compression  based  on  theoretical  analysis  and  1 D  calculations  as  described  in  the  proposal  dated  8- 
Aug-94. 


14.  SUBJECT  TERMS 

15.  NUMBER  OF  PAGES 

34 

Nil 

16.  PRICE  CODE 

N/A 

17.  SECURITY  CLASSIFICATION 

OF  REPORT 

18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

19,  SECURITY  CLASSIFICATION 

OF  ABSTRACT 

20.  LIMITATION  OF  ABSTRACT 

UNCLASSIFIED 

UNCLASSIFIED 

UNCLASSIFIED 

UL 

NSN  7540-01-280-5500 

Standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  Std.  239-18 

298-102 

INVESTIGATORS 


A.  M.  Buyko 
O.  M.  Burenkov 
F  Garanin 
- - V.  F.  Yermolovich 


V.  N.  Mokhov 


A.  I.  Startsev 


A.  I.  Sushko 
<^^V.  B. 


Yakubov 


Table  of  Contents. 


Table  of  Contents. . 3 

1.  Introduction . 4 

2.  Selection  of  physical  scheme . 5 

3.  Zero  -  dimensional  computation  setting  -  up  and  results . 6 

3.1  Computations  without  account  of  compressing  shell 

decebration . . 6 

3.2  Computations  with  account  of  compressing  shell 

deceleration . 12 

4.  Selection  of  working  material . 12 

5.  Results  of  one  -  dimensional  gas  -  dynamical  computations  and 

comparison  with  zero  -  mentional  computations . 20 

6.  Discussion  of  results . 31 

References . 34 


This  report  considers  capabilities  of  the  system  with  compression  of 
low  -  dense  plasma  heated  to  temperature  T0  =  0. 1^0.3  keV  with  a 
relatively  slow  liner  (V0  =  1-^3  cm/ps)  to  obtain  megajoule  amounts  of  X 
radiation. 

The  report  does  not  consider  a  specific  facility  design,  issues  of  plasma 
pre  -  heating  and  radiation  withdrawal.  At  this  stage  the  principal  objective 
is  to  find  out  possible  values  of  the  system  parameters  basing  on 
estimations,  analytic  and  one  -  dimensional  numerical  computations  which 
will  allow  to  consider  the  issue  whether  or  not  it  is  reasonable  to  develop 
this  direction  of  research. 
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2.  Selection  of  physical  scheme. 


Selection  of  the  physical  scheme  with  quasi  -  spherical  plasma  compression 
with  a  liner  was  discussed  in  the  first  part  of  the  report  [%].  It  is  assumed 
H  that  quasi  -  spherical  (not  necessarily  hydrogen)  plasma  compression  can  be 

made  using  a  system  similar  to  that  intended  for  reachment  of  DT  -  plasma 
ignition  threshold  at  plasma  chamber  compression  by  a  liner. 

We  will  consider  a  system  where  kinetical  energy  of  the  spherical  shell 
is  converted  to  X  radiation  energy  at  implosion  of  intermediate  “working 
substance”,  i.  e.  plasma,  heated  to  temperature  Tq  —  0.2  keV  at  the  initial 
spherical  shell  radius  of  Rq  =  6  cm.  Such  selection  of  initial  plasma 
parameters  prior  to  discussion  of  a  pre  -  heating  technique  is  arbitrary  to 
some  extent.  We  oriented  to  typical  plasma  parameters  in  the  MAGO 
system  chamber  which  can  be  considered  as  one  of  possible  pre  -  heating 
variants. 

In  our  computation  initial  kinetic  energy  of  the  shell  was  selected  equal 
to  Eq  =  20  MJ  which  corresponds  to  capabilities  of  appropriate  EMG. 

The  report  presents  similarity  parameters  allowing  to  update  basic 
system  characteristics  for  other  initial  conditions.  Initial  plasma  density  and 
shell  velocity  were  selected  so  that  maximum  plasma  temperature  be 
reached  at  compression  5  =  p/p0  ^  1000  which  seems  to  be  quite  real  for 

spherical  systems  since  shell  radius  will  vary  no  more  than  by  a  factor  of  10. 
Reachment  of  high  compressions  will  be  impeded  with  shell  motion 
instability.  The  possibility  to  obtain  X  radiation  with  quanta  energy  of 
hv~  5  4  10  keV  at  shell  velocity  of  Vq  ~  34-5  cm/ps  is  based  on  absence  of 
radiation  -  material  equilibrium  due  to  low  plasma  density.  For  low  -  z 
plasma  at  temperature  ~  5  4  10  keV  the  basic  radiation  mechanism  is 
Thus,  the  system  under  consideration  must  be 
spectrum  X  radiation  source  with  characteristic  emission 

«  0.2  ps. 

* 


bremsstrahlung. 

bremsstrahlung 

time  T  =  0.1 — - 

V0 


* 


3.  Zero  -  dimensional  computation  setting  -  up  and 

results. 


3.1  Computations  without  account  of  compressing  shell 

deceleration. 


To  find  out  principal  features  of  the  processes  under  consideration, 
treat  the  problem  of  adiabatic  compression  of  plasma  heated  to  initial 
temperature  Tq  by  a  shell  moving  at  constant  velocity  V0.  Plasma  is 
assumed  to  be  fully  ionized  and  low  -  dense,  i.e.  there  is  no  radiation  - 
material  equilibrium  owing  to  very  low  optical  thickness  of  plasma  under 
consideration.  Assume  that  plasma  energy  losses  to  radiation  are  described 
by  the  formula  for  bremsstrahlung  [2] 


q  =  1.97-105-p^-,/7;  \MJ/g^s} 
A 


(3.1) 


where  p  -  plasma  density  [g/cm3], 

Te  -  electron  temperature  [keV], 

Z  -  ion  charge, 

A  -  atomic  weight. 

Write  energy  balance  for  electron  and  ion  plasma  components  (per  lg 
of  material)  in  the  form: 


d&i  =  _p  dp^_  c  (Ti~Te) 

dt  dt  1  xie 

^L  =  -P  dp-‘  ,  c- ft*  -Tg)  a 
dt  dt  1  x  ie 

Here: 

=  CiTi  -  internal  energy  of  the  ion  component, 

Ce  =  ZC{Te  -  internal  energy  of  the  electron  component, 

Pj  =  —  C{pT{  -  ion  pressure, 

3 

2 

Pe  =  —ZCipTe  -  electron  pressure, 

Ci  -  plasma  ion  component  heat  capacity, 

Xia  —  electron  and  ion  temperature  smoothing  time, 
q  -  deceleration  energy  losses. 


(3.2) 


(3.3.) 
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If  temperature  smoothing  time  is  much  less  than  compressing 


characteristic  times 


Tie  «  X  = 


R, 

VJ 


then  the  adiabaticity  condition  is 


satisfied  for  the  equation  system  3.2,  3.3  which  can  be  written  in  the  form 

<?*ie 


T‘"Te>  +  l)c. 

In  this  case  the  equation  system  reduces  to  the  equation 


(3.4) 


d 

dt 


{  rp  \ 

vp^y 


P  '  Ci(z  +  1) 


(3.5) 


Introduce  the  following  notations: 
p0  -  initial  plasma  density, 

T0  -  initial  plasma  temperature, 

5  =  — —  (relative  density  or  compression  degree  of  plasma), 

Po 


S  =  e 


_  £f 

p2'3  T0 


-  a  dimensionless  value  whose  logarithm  determines 


electron  entropy, 

Rq  -  initial  density  radius, 

Vq  ~  shell  motion  velocity, 

_  rq 

X  = -  -  characteristic  compression  time. 

V0 

When  shell  moves  at  a  constant  velocity  the  compression  degree  is 
described  with  the  formula 

l 

5 3  =  ■—  (3.6) 

x  -  t 

Eq.  (3.5)  with  account  of  (3.1)  and  (3.6)  reduces  to  the  dimensionless  form 

dVs 


d5 


=  -aS-2/3 


(3.7) 


where 


a  = 


224.5R0p0z‘ 


v0Vvi(z  +  i) 

The  solution  of  Eq.  (3.7)  is  of  the  form 


(3.8) 


8 


s  =  (l  +  3af(l~- 


3  a 


5V3 


\2 


+  3d  J 

From  this  expression  it  is  seen  that  when  compression 


(3.9) 


5.  = 


1  +  3a 
3d  / 


(3.10) 


is  reached  plasma  temperature  vanishes.  The  corresponding  time  is 

5l/3  - 1 


t.  =  x 


5f 


(3.11) 


Using  formulas  (3.6),  (3.10),  (3.11)  one  can  represent  (3.9)  in  the  form 

2 


2/3 

f 

i  — 

f- 

.  1/3^ 

5f-l 

i  — 

V 

J 

(3.12  a) 


or 


s(t)  = 


-2  rt  t  2 


t: 


V  X  -  t 


(3.12  b) 


Plasma  density,  temperature  and  pressure  are  described  by  the  formulas 

_3 


P  -  Po5  -  Po 


(x-t)3' 


(3.13) 


T  =  T0S2/3s(8)  =  T„ 


(x-t)2 


(3.14) 


P  =  P„82/3s(6)=  P0 


(x-t)5 


<t). 


where  P0  =  —  (z  +  l)c<PoT0  -  initial  plasma  pressure.. 


Plasma  energy  is 


A  jr 

E  =  (z  +  l)CjT  — p0R03 


(3.15) 


or 
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* 


E  =  E052'3  ■  s(8)  =  Ec 


(t  - 1 ) 


<t) 


(3.16) 


where  £0  -  initial  plasma  energy. 

The  work  of  entire  plasma  compression  per  time  unit  is  of  the  form 


A  =  ^-5-s(8)=2£(  1 


<0 


T  (T-t) 

Intensity  of  plasma  radiation  from  the  whole  compressed  volume  is 


(3.17) 


•  2E0S4/3-sV2(5)  (T-t.)  t3 

e  =  T  ^3TT“  2E°-7T'(x-ty 


s1/2(t )  (3.18) 


The  integral  values  corresponding  to  (3.17)  and  (3.18)  (shell  work  on  plasma 
and  radiation  energy  released  from  plasma)  are  described  by  formulas 


A( 5)  =  E( 


2/3 


(5f  - 1) 


(82/3-1)-i^(5-1)^(54/3-l) 


2Sf 


(3.19) 


and 


Q(5)=EC 


2/3 


(8/3  -  lj 


38  f 


(3.20) 


Using  formulas  (3.15)  -  (3.18)  one  can  obtain  an  expression  for  maximum 
values  of  plasma  parameters:  *) 
maximum  temperature 


maxT  «  0.0625T0  -5f 


(3.21) 


is  reached  at  5  =  0.1258*; 
maximum  work  flow 


•  E 

max  A  «  0.691  —  5, 

T 


is  reached  at  8  =  0.2165*; 
maximum  radiation  intensity 


maxQ  «  0.164 


E 


0  5, 


(3.22) 


*) 


In  the  approximate  formulas  given  below  it  is  assumed  that 


6.'/J  »  /. 


is  reached  at  8  =  0.5125*. 

From  formulas  (3.19),  (3.20)  one  can  derive  expression  for  limiting 
values  of  work  done  on  plasma  and  emitted  energy(at  the  time  when 

5  =  5*) 


-A*  «  Q 


* 


—  5^3 
6 


(3.23) 


Figs.  3.1  and  3.2  give  plots  illustrating  formulas  (3.14)  -  (3.20).  In  the 
selected  scales  the  plots  are  of  universal  character  and  practically  do  not 
depend  on  particular  value  of  the  parameter  8*.  Temperature  was 
therewith  normalized  to  max  T ,  differential  works  and  radiations  to 

maxQ,  integral  works  and  radiation  flows  and  plasma  energy  to  A(5*). 


Basic  values  vs  time  when  shell  moves  at 
constant  velocity. 


The  principal  feature  of  the  obtained  solution  is  that  beginning  from 
certain  time  plasma  radiation  intensity  becomes  higher  than  work  flow  at 
the  boundary  and  plasma  begins  to  cool  dramatically.  In  so  doing, 
irrespective  of  infinite  shell  energy  (absence  of  deceleration)  finite  work 
(3.23)  is  being  done  on  plasma.  The  X  radiation  pulse  has  a  stretched  leading 
and  very  steep  back  front  (Fig.  3.2).  Maximum  radiation  intensity  is  reached 
when  plasma  compression  is  4  times  higher  than  at  maximum  temperature, 
plasma  temperature  being  reduced  down  to  T  =  0.4  max  T . 

Consider  some  quantitative  estimations.  Assuming  that  plasma 
compression  by  the  shell  is  constant  up  to  8  =  1000,  select  8*  =  8000, 

then,  according  to  (3.21),  at  8  =  1000  maximum  temperature  T  25T0 , 


internal  plasma  energy  E  «  25 EQ ,  differential  work  flow 


A  =  500 


•  E0 

radiation  intensity  Q  =  500  — are  reached.  By  that  time  plasma  emits 

X 

energy  Q  =  21E0,  while  the  ultimate  value  of  emitted  energy,  according  to 

(3.23),  is  Q*  =  67.E0.  Maximum  work  flow  is  reached  at  8  =  1720  and  is 

•  Eq  •  E0 

A  =  550  — ,  at  that  time  radiation  intensity  is  Q  =  830 —  at 

T  T 

temperature  T  =  23  T0.  By  that  time  energy  Q  =  32E0  is  emitted. 


E 


Maximum  radiation  intensity  Q  =  1300--^  is  reached  at  compression 

X 

8  =  4000,  at  temperature  T  =  10T0.  By  that  time  energy  Q  =  55£0is 
emitted. 


3.2  Computations  with  account  of  compressing  shell 

deceleration. 


To  account  shell  acceleration,  complete  Eq.  (3.5)  with  the  equation  of 
shell  motion 

dV 


m 


dt 


=  -4t zRzP, 


(3.24) 


where  m  -  shell  mass,  V  -  shell  velocity.  Introduce  dimensionless  velocity 

Al 
wn 


V  ,  r 

U  — -  and  two  dimensionless  parameters  %  =  1  —  VS  and  k  = 


v ; 


where  W0  = 


mV' 


-  initial  shell  energy.  Transform  Eq.  (3.5)  to  the  form 


dt 


d%  3  a82/3 

Eq.  (3.24)  with  account  of  (3.25)  is  reduced  to  the  form 


(3.25) 


du 

d% 


A 

3a 


<5 1/3  (l  -  if 


(3.26) 
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Formula  (3.6)  has  to  be  replaced  with  the  equation 

dS1/3  u 


(3.27) 


dx  3  a 

Equation  system  (3.26),  (3.27)  yields  the  parametric  solution  to  the  problem 
and  Eq.  (3.25)  sets  correlation  of  the  parameter  %  to  real  time.  This  equation 

system  is  equivalent  to  the  second  -  order  linear  equation 

whence  it  is  seen,  that  another  dimensionless  combination  of  initial 
parameters  which  determines  the  course  of  the  solution  is  the  expression 


Ke  = 


Er 


E 


Wn  •  9a2  Wn 


°--5f  «6 


& 

Wn 


(3.28) 


Fig.  (3.3)  gives  results  of  numerical  solution  to  Eqs.  (3.26),  (3.27)  for  a 
set  of  Ke  values  at  8*  =  8000.  At  Ke  =  12  (W0  =  0.5Q*  complete  shell 
deceleration  takes  place  at  %  =  1.  Apparently,  in  this  case  the  maximum 
system  efficiency  factor  is  realized:  entire  shell  energy  is  converted  to 
radiation.  At  lower  values  of  Ke  plasma  does  not  have  time  to  extract 
whole  shell  energy  due  to  rapid  cooling,  at  high  values  of  Ke  plasma  does 
not  have  time  to  convert  kinetic  energy  of  the  shell  to  radiation.  It  should 


be  noted,  that  the  parameters  Ke  ~ 


Ml 

mp0 


does  not  explicitly  depend  on 


initial  shell  velocity. 

Cite  the  formulas  for  computation  of  physical  parameters  through 
solution  of  Eqs.  (3.26,  3.27): 

E  =  E051/3(l-x)2, 

T  =  T082/3(1  -x)2, 


Q  = 


A  = 


2E„84/3(1-x) 
(6f-l)  ’ 


5(1  -  if  v- 


X 

2  E, 


(3.29) 


A  =  W0(l-u2), 

Q  =  A  +  E0  -  E. 

Figs.  3.4  -  3.7  give  plots  of  these  values  computed  for  several  values  of 
Ke,  with  the  parameter  a  being  selected  so  that  maximum  temperature  be 
reached  at  5  =  1000.  In  these  plots  temperature  is  given  in  units  T0, 

E0 

energy  and  work  in  units  E0 ,  differential  flows  in  units  — . 


V/Vo 


16 

Solution  of  equation  system  (3.  26),  (3.  27) 
for  Ke=3,  6*  =121  67 


4 


* 


* 


5 


■c 

■or 


■< 

■Q? 


* 
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Solution  of  equation  system  (3.  26),  (3.  27) 
for  Ke=6,  6*  =17985. 


t/tau 
Fig.  3.  6 
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4.  Selection  of  working  material. 


* 


* 


* 


The  analysis  performed  in  the  previous  section  imposes  no  explicit 
limitations  on  the  value  z.  The  limitation  on  z  follows  from  the  applicability 
condition  of  formula  (3.1)  which  describes  plasma  radiation  within  the 
temperature  range  providing  not  only  total  material  ionization,  but  also 
small  contribution  of  bound  -  free  transition  radiation.  For  the  initial 

temperature  T0  =  0.2  keV  only  hydrogen  and  helium  satisfy  these 
requirements.  Description  of  high  Z  materials  requires  specification  of  the 
above  model.  We  did  it  using  the  interpolation  formulas  given  in  paper  [3]. 

With  account  of  energy  expenditure  on  ionization  the  equation  of  state 
of  plasma  was  written  in  the  form 


r  =  |  (z  + 1)  •  c  <  •  P  •  t 

8  =  C;  •  (z  +  l)-  T 
=  cj-(z  +  l)+ej-T.^  + 


dQz 

dT 


(4.1) 


where  ct  -  heat  capacity  of  ion  ideal  gas, 
z(T)  -  mean  ion  charge, 

QZ(T )  -  expenditure  of  energy  on  ionization. 

In  the  approximation  accepted  in  [3]  mean  ion  charge  depends  only  on 
temperature,  in  connection  with  this  internal  energy  and  heat  capacity  of 
plasma  in  (4.1)  are  functions  of  temperature  only.  For  simplicity  the 

dependency  QZ(T )  was  taken  in  the  form  QZ(T)  =a(z(T^3 ,  with  the 

coefficient  a  being  chosen  so  that  proper  energy  of  total  atom  ionization  be 
secured. 

With  account  of  plasma  heat  capacity  dependence  on  temperature  Eq. 
(3.5)  takes  the  form 

d  f  T  ' 

dt  v  p2/3  y 


f  m  dz  dQ 

C:T  -  +  - - 

2  V  dT  dT 


vp^3/ 


dp 
pdt  c  p 


2/3 


(4.2) 


The  first  addend  in  the  right  -  hand  side  accounts  consumption  of  energy 
on  ionization,  the  second  -  energy  losses  on  radiation.  Instead  of  formula 
(3.1)  the  following  expression  was  used  for  q 

q  =  Pz(T)-pz(T),  (4.3) 

where  z  (T)  and  Pz  (T)  are  described  by  the  interpolation  formulas  from  [3] 
Eq.  (4.2)  was  numerically  solved  simultaneously  with  equations  of  shell 
motion  (3.24).  The  initial  parameters  of  the  system  were  selected  so  that 
maximum  plasma  temperature  be  reached  at  compression  8  ~  1000  and 
whole  plasma  energy  be  emitted  by  the  time  of  complete  shell  deceleration 
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(the  mode  with  maximum  efficiency).  The  initial  shell  radius  was  taken 
equal  to  R0  =  6  cm,  the  initial  shell  energy  -  20  W0  »  MJ.  The  computed 
results  are  given  in  Table  4.1  and  in  plots  of  Figs.  4.1  through  4.8. 

The  computations  show  that  at  the  initial  temperature  T0  «  0.2  keV 
acceptable  results  are  obtained  only  for  low  z  (He,  Li).  As  z  increases,  with 
the  initial  temperature  being  fixed,  shell  velocity  needed  to  implemerit  the 
selected  mode  increases  and,  accordingly,  shell  mass  decreases.  Already  for 
nitrogen  (z  =  7)  at  T0  =  0.2  keV  it  is  impossible  to  realize  the  mode  with 

complete  shell  deceleration  at  reasonable  values  of  Vq  (see  computations 
Nos.  3,4  in  Table  4.1).  If  maximum  temperature  is  therewith  reached  at 
5  ~  1000,  then  the  mode  with  low  efficiency  (~  10  %)  is  realized  and 
maximum  plasma  temperature  is  too  low  (computation  No.  3,  Table  4.1).  The 
complete  shell  deceleration  mode  is  realized  at  plasma  compression 

5  ~  2  •  104  (computation  No.  4,  Table  4.1).  It  should  be  noted  that  transition 
from  mode  to  mode  occurs  at  relatively  small  change  in  initial  density. 

Computations  Nos.  7,  8  of  Table  4.1  show  that  for  materials  with 
z  =  7...  10  initial  plasma  temperature  should  be  Tq~  0.5. ..0.9  keV. 

Mention  one  more  feature  of  the  solution  using  formula  (4.2)  for 
nitrogen  and  neon.  At  plasma  cooling  down  to  temperature  ~1  keV  on  the 
• 

dependency  Q(t )  the  second  acute  maximum  appears  which  is  related  to 

abrupt  increase  in  energy  losses  to  radiation  within  this  temperature  range, 
but  energy  emitted  in  this  peak  is  negligible. 

Keeping  in  mind  the  possibility  to  use  the  MAGO  chamber  for  plasma 
preheating,  we  computed  the  system  with  the  central  electrode  of  0.1  R0 
(computations  Nos.  9,  10  in  Table  4.1).  Availability  of  the  electrode  led  to  a 
small  decrease  in  maximum  temperature  as  compared  to  a  similar 
computation  of  the  system  with  free  center  (computation  No.  1). 
Computation  No.  11  was  performed  for' comparison  with  computation  No.  7 
of  Table  5.1. 

The  computed  results  show  that  for  “working  substance”  it  is  most  real 
to  use  low  z  material:  helium  or  even  hydrogen.  High  z  materials  appear 
more  preferable  from  the  source  intensity  increase  standpoint,  but  impose 
more  stringent  requirements  on  the  initial  shell  velocity  and  the  initial 
plasma  temperature. 
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5.  Results  of  one  -  dimensional  gas  -  dynamical 
computations  and  comparison  with  zero  - 
dimentional  computations. 


The  afore  -  stated  results  allow  to  estimate  the  values  of  the  initial 


system  parameters  close  to  optimal.  The  search  for  these  values  by  the 
results  of  direct  numerical  experiment  seems  to  be  if  not  hopeless  then 
extraordinarily  expensive  and  tiresome  affair.  Moreover,  as  the  comparison 
with  the  results  of  one  -  dimensional  computations  shows,  the  used  method 
exhibits  acceptable  accuracy  for  estimation  of  expected  system 
characteristics. 

Since  at  this  phase  of  the  works  selection  of  specific  design  is  not 
suggested  then  the  one  -  dimensional  computations  performed  are  of  model 
character  and  must  be  considered  as  a  further  justification  and  specification 
of  the  above  results.  In  addition  to  more  accurate  computation  of  plasma 
implosion  dynamics,  the  one  -  dimensional  computations  accounted  effect  of 
the  processes,  such  as  electron  and  ion  heat  conduction  and  plasma 
magnetization. 

System  No.l  from  Table  4.1,  i.e.  helium  plasma  with  initial  density 
p0=  1.3  g/cm3  and  initial  temperature  T0  =  0.2  keV  in  spherical  volume  of 
jR0  =  6  cm  initial  radius,  was  used  as  the  basis  for  the  numerical 
experiments.  Plasma  was  compressed  with  a  shell  m  =  45  g  in  mass  at  initial 
velocity  V0  =  3  cm/ gs.  When  modeling  the  state  of  the  electrically 
“exploded”  liner  we  set  the  initial  shell  density  ~  0.1  g/cm3.  Plasma 
magnetization  was  modeled  by  “frozen  -  in”  magnetic  field  with  the  initial 


ratio 


h; 


8nPn 


=  0.3 .  The  computed  results  are  given  in  Table  5.1. 


Computation  No.  1  shows  that  the  model  considered  in  Section  4  is 
sufficiently  accurate  to  perform  preliminary  estimations  of  the  system 
parameters.  It  represents  maximum  plasma  temperature  well.  Reduction  in 
the  fraction  of  shell  energy  converted  to  radiation  energy  down  to  78  % 
(instead  of  100  %)  is  related  to  account  of  expenditure  of  energy  on  shell 
compression.  This  energy  fraction  is  mainly  spent  for  implosion  of  the  shell 
portion  “evaporated”  by  radiation  whose  mass  by  the  end  of  the  process  is 
considerably  higher  than  plasma  mass  and  essentially  depends  on  shell 
material. 

Maximum  plasma  compression  and  maximum  radiation  intensity 
reduce  for  the  same  reason.  The  model  computations  of  Section  4 
overestimate  the  maximum  radiation  intensity  value  by  a  factor  of  ~  2  . 

Plasma  magnetization  (computation  No.  2)  weakly  affects  maximum 
temperature,  but  somewhat  reduces  compression  and  the  fraction  of  energy 
converted  to  radiation. 

Account  of  electron  and  ion  heat  conduction  is  represented  by 
computations  Nos.  3,  4.  As  compared  to  computation  No.  2,  they  show  that 


“freezing”  of  electron  heat  conduction  with  magnetic  field  is  needed  while  of 
ion  heat  conduction  desirable. 

Computations  Nos.  5,6  accounted  the  effect  of  the  central  electrode,  it 
affected  maximum  plasma  compression  and  radiation  output  but  not 
maximum  temperature. 
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6.  Discussion  of  results. 

The  computations  made  show  basic  possibility  to  obtain  a 
bremsstrahlung  spectrum  X  radiation  source  of  characteristic  temperature 
5  v  7  keV,  emission  time  x  «  0.2  ps  and  energy  conversion  factor  >  50  % 
through  spherical  implosion  of  low  -  density  helium  plasma  preheated  up  to 
temperature  Tq  =  0.2  keV 

For  practical  employment  of  such  a  source  it  is  necessary  to  solve  the 
problem  of  radiation  extraction  from  cavity.  It  is  easy  to  verify  that 
radiation  extraction  directly  through  the  shell  is  impossible.  Indeed,  “optical” 
thickness  of  the  shell  in  the  final  state  for  the  system  considered  in  par.  5  is 
pAR  ~10  g/cm2  even  if  radiation  range  in  the  shell  is  purely  Compton 

2 

(l  *  5/p  cm),  then  radiation  will  be  already  attenuated  by  a  factor  of  e  . 

Moreover,  as  the  computations  show,  the  use  of  low  z  materials  in  the 
shell  noticeably  reduces  (by  ~  30  %)  the  shell  energy  fracture  converted  to 
radiation. 

It  is  believed  that  radiation  can  only  be  extracted  through  Be  window 
of  ~  0.1  R0  radius  *)  located  on  the  system  axis.  “Optical  thickness”  of  such  a 
window  will  be  ~  0.5  -r  1  g/cvn2  with  mass  0.5  t  1  g.  Radiation  absorbed  by 
the  high  -  z  shell  is  converted  to  equilibrium  X  radiation  whose 
temperature  reaches  ~  0.3  keV  in  the  end  of  the  process.  This  energy  can 
heat  up  the  window  to  the  “transparency”  level  to  secure  fairly  efficient 
extraction  of  hard  radiation.  Apparently,  it  is  this  aspect  of  the  problem 
that  will  eventually  determine  needed  energetics  and  optimal  parameters  of 
the  system  .  Discussion  of  this  problem  in  more  detail  is  beyond  the  scope  of 
the  paper. 


*)  For  the  system  with  the  central  electrode  it  will  be  an  annular  window 
approximately  of  the  same  area. 
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